Spin Orientation of Holes in Quantum Wells 
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This article reviews the spin orientation of spin-3/2 holes in quantum wells. We discuss the Zee- 
man and Rashba spin splitting in hole systems that are qualitatively different from their counterparts 
in electron systems. We show how a systematic understanding of the unusual spin-dependent phe- 
nomena in hole systems can be gained using a multipole expansion of the spin density matrix. As 
an example we discuss spin precession in hole systems that can give rise to an alternating spin 
polarization. Finally we discuss the qualitatively different regimes of hole spin polarization decay 
in clean and dirty samples. 



I. INTRODUCTION 

Spin electronics is a quickly developing research area 
that has yielded considerable new physics and the 
promise of novel applications The electrons in the 
conduction band of common semiconductors like GaAs 
are characterized by a spin-1/2. Holes in the topmost 
valence band, on the other hand, have an effective spin- 
3/2 (Ref. 0) which gives rise to many novel features 
that are not present in the conceptually simpler case of 
spin-1/2 electron systems. 

In this article we review some of the intriguing phe- 
nomena related with the spin orientation in spin-3/2 hole 
systems. We begin in Sec. [TT] with a brief review of the 
Luttinger Hamiltonian which forms the foundation for a 
theoretical description of spin-3/2 hole systems in cu- 
bic semiconductors like GaAs. In Sec. [Hi] we discuss 
the anisotropic Zeeman splitting of two-dimensional (2D) 
hole systems. While a magnetic field B perpendicular 
to the 2D plane gives rise to a large Zeeman splitting, 
the splitting in an in-plane magnetic field is greatly sup- 
pressed. However, it can be tuned, e.g., by varying the 
thickness of the quasi-2D system by means of external 
gates. At B = 0, the Rashba spin-orbit coupling in an 
inversion-asymmetric 2D system is characterized by an 
effective magnetic field oriented in the 2D plane. As dis- 
cussed in Sec. IIV1 the resulting spin splitting in 2D hole 
systems behaves thus very similar to the Zeeman splitting 
in an external magnetic field. In particular, Rashba spin 
splitting can likewise be tuned by varying the thickness 
of the quasi-2D system. 

In Sec. [V] we review the multipole expansion of the 
spin density matrix that provides a more systematic un- 
derstanding of the unusual spin-dependent phenomena 
in hole systems. As a first application of this general 
approach, we discuss in Sec. IVII the multipole moments 
induced in a 2D hole system by an in-plane magnetic 
field. Next we use the multipole expansion to discuss in 



Sec. IVIII the spin precession in 2D hole systems which 
turns out to be qualitatively different from the more fa- 
miliar case of spin precession in spin-3/2 hole systems. 
For example, the hole spin polarization and the higher- 
order multipoles can precess due to the spin-orbit cou- 
pling in the valence band, yet in the absence of exter- 
nal or effective magnetic fields. Finally, we discuss in 
Sec. I Villi the spin polarization decay in hole systems. 
Here, an important parameter is the product of the pre- 
cession frequency il times the momentum relaxation time 
t p . Qualitatively different regimes can be distinguished 
for dirty samples with f2r p <^ 1, weak-scattering samples 
with Qt p 3> 1 and ballistic systems with r p — > oo. In 
Sec. HXI we summarize our results. 



II. 



SPIN-3/2 HOLE SYSTEMS: THE 
LUTTINGER HAMILTONIAN 
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In a tight-binding picture, the electrons in the con- 
duction band of common semiconductors like GaAs are 
described by s-like atomic orbitals 0,01, see Fig. Q] (left). 
Taking into account spin, the electrons have a total an- 
gular momentum j = 1/2 that behaves analogously to a 
spin s = 1/2. Holes, on the other hand, are described by 
p-like atomic orbitals. Taking into account spin, the holes 
have a total angular momentum j = 3/2 and j = 1/2. 
The atomic spin-orbit coupling separates the j = 3/2 
from the j = 1/2 states by a spin-orbit gap Ao, see 
Fig. [1] (center) , so that we can associate an effective spin 
j = 3/2 with the states in the topmost valence band. 
At nonzero wave vectors k the fourfold degenerate states 
j = 3/2 split further into so-called heavy-hole (HH) and 
light-hole (LH) states. If we choose the quantization axis 
of angular momentum j parallel to k, the HH (LH) states 
corresponds to m — ±3/2 (m = ±1/2). (We ignore here 
small fc-linear terms that couple HH and LH states Q.) 

In the following we consider quasi two-dimensional 
(quasi-2D) systems in the xy plane. Subband quanti- 
zation in quasi-2D systems correspond to standing waves 
in the z direction so that we get a splitting of HH and LH 
states called the HH-LH splitting, even for in-plane wave 
vector fen = 0, see Fig. [TJ (right). For 2D systems grown 
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FIG. 1: Qualitative sketch of the band structure of GaAs close 
to the fundamental gap. (Left) The electrons in the conduc- 
tion band are described by s-like atomic orbitals whereas the 
holes in the the topmost valence band are described by p- 
like atomic orbitals. (Center) Taking into account spin, the 
conduction band states are characterized by a total angu- 
lar moment j = 1/2, whereas in the valence band we have 
j — 3/2 and j = 1/2 states seperated by a spin-orbit gap 
Ao. (Right) In quasi-2D systems the j = 3/2 states split into 
HH and LH states. For electron and LH states the spin is 
essentially decoupled from the orbital motion, whereas HH 
states are predominantly made of m — ±3/2 states with a 
spin quantization axis perpendicular to the 2D plane. 



on a high-symmetry (001) or (111) surface, the k\\ = 
eigenstates are pure HH or LH eigenstates (if we neglect 
the small A:-linear terms). For any other surface, e.g., the 
experimentally important (113) surface we obtain 

a weak coupling of HH and LH states even at kxx = 
which is caused by the terms of cubic symmetry in the 
Hamiltonian (see below). 

For k\\ ^ we get a mixing of HH and LH states. 
However, for typical sample parameters (well width and 
Fermi wave vector /of), the wave vector k z characteriz- 
ing the standing wave perpendicular to the 2D plane is 
generally much larger than kp, k z 3> kp, so that we can 
interpret the subband states even for fcy ^ as HH or LH- 
like. Often, only the lowest HH subband is occupied, i.e., 
the occupied subband states are predominantly made of 
to = ±3/2 states assuming a spin quantization axis per- 
pendicular to the 2D plane. It is this fact that lies at 
the heart of several anomalous properties of 2D hole sys- 
tems that we are going to discuss in subsequent sections. 
Of course, at magnetic field B = the m = ±3/2 and 
m = —3/2 states contribute with equal weight, i.e., as ex- 
pected for a paramagnetic material, the effects discussed 
here do not give rise to a permanent magnetic moment. 

For a more explicit and detailed discussion of spin ori- 
entation we need to use the Hamiltonian appropriate for 
these hole systems. Generally, the j = 3/2 hole states 
in the topmost valence band of cubic semiconductors like 
GaAs are described by the Luttinger Hamiltonian [|J 
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the Luttinger parameters, J, are the 4x4 matrices for 
angular momentum j = 3/2 (see, e.g., Ref. [H), we have 
{A, B} — 7}(AB ± BA) and cp denotes cyclic permuta- 
tion. In Eq. |T]) we neglected the terms linear and cubic 
in k which contribute to the spin splitting due to bulk 
inversion asymmetry [4]. We remark that for the effects 
discussed here these terms are of minor importance. Us- 
ing an explicit matrix notation, Ti.^ can be written in the 
form 
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Here we have expressed J; and TLl in a basis of j = 
3/2 angular momentum eigenfunctions in the order to = 
±3/2, ±1/2, -1/2, and -3/2. 

The notation used in Eq. {l} reflects the cubic symme- 
try of the crystal structure. An alternative formulation 
of the Luttinger Hamiltonian was proposed by Lipari and 
Baldereschi Q 
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where 7 = (272 ± 373)/5 and J = {J x ,Jy 
first two terms in Eq. ([4]) have spherical symmetry with 
the first term being diagonal in spin while the second 
term can be interpreted as a spherical spin-orbit cou- 
pling within the j — 3/2 space, see also Eq. (TT9l) below. 
Finally, He. represents the anisotropic terms with cubic 
symmetry [8] which will be given in Eq. (|2"D|) below. Usu- 
ally the terms in 7i c are small. Neglecting these terms 
corresponds to the spherical approximation H. s of TLl- 
Using for TL S an explicit matrix notation as in Eq. we 
get P s — P and 
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where k± = k x ± ik y . In the spherical approximation, 
the energy dispersions for the HH and LH states are 



where toq is the mass of free electrons, 71, 72, and 73 are 
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For quasi-2D systems it is often advantageous to use 
an alternative decomposition Hl = H ax +H' c , where H ax 
has axial symmetry with the symmetry axis n chosen per- 
pendicular to the 2D plane. Neglecting H' c corresponds 
to the axial approximation [J, |9j. For n parallel to [001] 
we obtain for H ax , using an explicit matrix notation as 
in Eq. 0, 



P n 
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Explicit expressions for other crystallographic orienta- 
tions of n are given in Refs. 0, H . For quasi-2D sys- 
tems the advantage of the axial approximation over the 
spherical approxiamtion lies in the fact that 7i ax captures 
the most important physics of different crystallographic 
directions n while both H s and H ax yield a rotational 
symmetry with respect to the axis n. 

We can readily see from Eq. ((2|) that Hl becomes di- 
agonal for k x = k y = 0, i.e., subband states at fe|| = 
are either pure HH (m = ±3/2) or pure LH (m = ±1/2) 
states. As discussed above, in a coordinate frame where 
the z axis points in a crystallographic direction other 
than the high-symmetry directions [001] or [111], Hl con- 
tains off-diagonal terms proportional to k z so that even 
at the subband edge fey = we get a mixing of HH and 
LH states. 



III. ANISOTROPIC ZEEMAN SPLITTING 

For spin-1/2 electron systems it is well-known that the 
Zeeman splitting, i.e., the response of the electron's spin 
to a magnetic field B, is usually essentially independent 
of the orientation of B [lG(. For electrons, the differences 
between the in-plane and perpendicular effective Lande 
factors <7j| and g^_ are most pronounced in narrow GaAs- 
AlGaAs quantum wells, where gjj and g* ± change sign and 
thus cross zero as a function of well width, as discussed in 
Refs. [0, [3- The in-plane anisotropy of the in-plane 
in low-symmetry geometries was discussed in Ref. [l3j . 
see also Ref. 

The situation in HH systems is qualitatively different. 
As suggested by Fig. [1] the response of a HH system to 
an in-plane magnetic field B\\ is suppressed because the 
effect of Bn competes with the rather rigid spin orienta- 
tion induced by the HH-LH splitting. A perpendicular 
magnetic field B±, on the other hand, is compatible with 
the B = spin orientation of the m — ±3/2 HH states so 
that we can have a large Zeeman splitting. The resulting 
anisotropy of g* was first discussed in Ref. [l4T ]. 

A more quantitative discussion needs to be based on 
the Luttinger Hamiltonian Hl with the Zeeman term 
Hz = — 2/iB kB • J added to it. Here k is the isotropic g 
factor. We neglect the anisotropic Zeeman term because 
its prefactor (often denoted q [2J) is usually much smaller 
than k. Using an explicit matrix notation, the Zeeman 



term reads 
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where B± = B x ±iB y . We can immediately read off from 
Eq. © that a perpendicular magnetic field B z gives rise 
to a Zeeman splitting A_Ehh = 6/ibkB z of the HH states 
whereas for LH states we get Ai?HH = 1\i^kB z . Also, 
we see from Eq. ([8]) that, in the presence of an in-plane 
magnetic field Bp = (B x , B y ,0), the Zeeman term Hz 
couples the two LH states with AE'lh = 4/iB/C-Bn, and it 
couples the HH states to the LH states [la ]. But there 
is no direct coupling between the HH states proportional 
to k, so that the Zeeman splitting of HH states in an 
in-plane magnetic field is suppressed [141 ]. 

A quantitative analysis based on perturbation theory 
shows [3] that, for fey = and neglecting the cubic terms 
H c , there is no Zeeman splitting of HH states linear in 
Bn , but the lowest-order splitting is proportional to B?> 
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E\ - E\, with E\ and E[ the energy of 
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obtain for fen > 
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As expected from Fig. [T] the energy denominators in Eq. 
([9|) reflect the competition between the spin orientation 
perpendicular to the 2D plane induced by the subband 
confinement and the (generally weaker) effect of the in- 
plane magnetic field B\\ that tends to orient the spins 
in-plane. The fen dependence in Eq. (|9b[) originates in 
the HH-LH coupling at k\\ > 0. It follows from Eq. © 
that we can tune the Zeeman splitting of HH systems in 
an in-plane magnetic field if we change the confinement 
in z direction by means of front and/or back gates. Al- 
ternatively, one can also change the HH-LH splitting by 
means of strain. 

When taking into account the cubic terms H c in the 
Luttinger Hamiltonian, the Zeeman splitting linear in Bv, 
(fe|| = 0) still vanishes for the high-symmetry surfaces 
(001) and (111). For other surfaces, we get a JSy-linear 
splitting the magnitude of which depends on the in-plane 
orientation of By relative to the crystal axes [liij . For 
an infinitely deep rectangular QW grown on an (mmn) 
surface we obtain in second order perturbation theory 
the following Zeeman terms acting within the space of 
the topmost HH subband 
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FIG. 2: Anisotropic effective <? factor g* of the hi subband 
for a GaAs-AIGaAs QW as a function of 8, t he ang le between 
[001] and the growth direction. Here x || [nn(2m)\, y || [110] 
and z II \mmn]. 



where 



9 XZ 
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Here the x direction corresponds to [nn(2m)], y corre- 
sponds to [TlO] and is the angle between [mmn] and 
[001], i.e. 9 = arccos(n/v / 2m 2 + n 2 ). We obtain for the 
coefficient /C 



/C = 



k(73 - 72) 



4 [(1 - C) 72 + C73] 



where 



C = sin 2 9 [3 



I sin 2 
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Note that these expressions are independent of the width 
of the QW. The importance of Eq. (fT0|) lies in the fact 
that the topmost subband in an (unstrained) QW is an 
HH subband so that often only this subband is occupied. 
In Fig. [2] we show g* x = g* y and g* z for the topmost 
HH subband in a GaAs-AIGaAs QW as a function of 
the angle 9. Figure [5] demonstrates that g* can be very 
anisotropic. For example, for the growth direction [113], 
g* z is more than a factor of four larger than g* x — g* y . 
For comparison, we remark that for the GaAs system 
considered in Fig. [2] we have g* zz = 6k w 7.2. Equations 
(fTT)|) and (fTTj) are applicable to a wide range of cubic 
semiconductors with results qualitatively very similar to 
Fig.d 

The anisotropic Zeeman splitting (JTUJ) has been studied 
experimentally by measuring the magnetoresistance of a 
high-mobility 2D hole system as a function of in-plane 
magnetic field -By [lg]. The sample was a 200 A wide 
Si-modulation doped GaAs QW grown on (113)A GaAs 
substrate. The left two panels of Fig.[3]show the resistiv- 
ity p measured as a function of B» for different directions 
of B|| and current I and for three different densities. For 
easier comparison we have plotted the fractional change 
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FIG. 3: Left and central panels: Fractional change in resistiv- 
ity p(_B||)/p(0) due to an in-plane B, measured at T — 0.3 K 
in a GaAs 2D hole system grown on a (113) substrate, for 
different directions of By and I and different 2D densities as 
indicated. The arrows mark B* as defined in the text. Right 
panel: Calculated density p+ in the upper spin subband as a 
function of Bu . 



p(i?ll)/p(0). It can be seen that log(p) shows a change in 
slope at a value of _By we call B* . In Fig. [3]£?* is marked 
by arrows. This magnetoresistance feature is related to 
a spin-subband depopulation and the resulting changes 
in subband mobility and intersubband scattering as Bn 
is increased [ij], Ha]- ^ ^ s remarkable that B* for the 
B || [332] traces is several Tesla smaller than for the [110] 
traces, but it is independent of the direction of I. This 
is strong evidence for the anisotropy of the in-plane g* . 
The experimentally observed anisotropy is qualitatively 
consistent with our self-consistently calculated results for 
the density p + of the upper spin subband as a function 
of -By , shown in the right panel of Fig. [3l The density p+ 
decreases much faster for B || [332] than for B || [1 10], in 
agreement with Fig. [21 The findings have been confirmed 
by Shubnikov-de Haas measurements probing directly the 
depopulation of the minority spin subband for different 
orientations of By [l9j . 

In Fig. [3] the measured B* is significantly smaller than 
the calculated By for a complete depopulation of the up- 
per spin-subband. We note that for our low-density sam- 
ples it can be expected that g* is enhanced due to the ex- 
change interaction and the spin polarization caused by -By 
[ToL [I?], H3] • These many-particle effects were not taken 
into account in our self-consistent calculations. However, 
they do not qualitatively affect the anisotropy of g* [l(| • 
See Ref. [2l| for a more detailed discussion of exchange- 
correlation effects in low-density hole systems on a (113) 
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IV. RASHBA SPIN SPLITTING IN 2D HOLE 
SYSTEMS 

If in a solid the spatial inversion symmetry is broken, 
we can have a spin splitting of the electron and hole states 
even in the absence of a magnetic field B. In quasi-2D 
semiconductor structures, the bulk inversion asymmetry 
(BIA) of the underlying crystal structure (e.g., a zinc 
blende structure) , and the structure inversion asymmetry 
(SIA) of the confining potential are usually the dominant 
contributions to the B = spin splitting (22J. While 
BIA is fixed, the so-called Rashba spin splitting [23| due 
to SIA can be tuned by means of external gates that 
change the electric field £ in the sample 0, [2f| . 

Here we want to focus on Rashba spin splitting. In 
quasi- 2D electron systems it is described by the term 



H C R = a£ z i (fc_CT+ — fe+<x_) = a£ z 




-ik-i 



ik- 




(12) 

where a± = \{px ± i& y ), a is a material-dependent pref- 
actor and £ z is an (effective) electric field that character- 
izes the broken inversion symmetry of the sample. Often 
it is instructive to write the Rashba term (I12|) in the form 
of a Zccman term 126] 



±B(k|, 



(13) 



where B(k\\) — 2a£ z (k y ,—k x ,0) is an effective magnetic 
field. It follows from the defintion of B(ky) that B ori- 
ents the electron spin in the plane of the 2D structure, 
perpendicular to the wave vector k|| . 

For j = 3/2 hole systems we need to replace the 2x2 
Pauli matrices a x and a y in Eq. (Ti"2"j) or (fT3"j) by the corre- 
sponding 4x4 matrices J x and J y for angular momentum 
j = 3/2, while the definition of the effective magnetic 
field S(kn) is valid also for hole systems, i.e., exactly as 
in electron systems, the field S(kii) tends to orient the 
hole spins in-plane. Here, many arguments on Zeeman 
splitting of hole states in an in-plane magnetic field carry 
over to the splitting in the effective field B(ky). In par- 
ticular, similar to Eq. ©, Rashba splitting of hole sys- 
tems is proportional to fcjj which is qualitatively different 
from the kn -linear Rashba spin splitting of electron sys- 
tems. Within the two-dimensional subspace of one HH 
subband, the Rashba term describing the cubic splitting 
is of the form 



Ttft = (3 h £ z i {k%a_ - k 3 _a^ 



(14) 



In electron systems, the coefficient a depends only on the 
material, but it is essentially independent of the geometry 
of the quasi-2D system [J]. In HH systems, on the other 
hand, we get in third-order perturbation theory for the 
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where a = 64/(97r 2 ) for an infinitely deep, rectangular 
QW, i.e., as expected from Fig. [TJ the prefactor (3^ de- 
pends on the HH-LH splitting which, in turn, depends 
on the geometry of the quasi-2D system. The detailed 
comparison of Eqs. © and (Ti"5)) reveals a subtle differ- 
ence between Zeeman splitting and Rashba spin splitting 
in HH systems. While Eq. ([9]) can be derived in second- 
order perturbation theory, each term in Eq. (|15j) contains 
two energy denominators 1/A. This is due to the fact 
that the effective field B depends on the wave vector kn 
whereas B is independent of kn . 

The functional form of the HH Rashba coefficient j3 h 
gives rise to a remarkable difference between Rashba spin 
splitting in electron and hole systems. It is well-known 
that Rashba spin splitting in electron systems is roughly 
linearly proportional to the electric field £ z that charac- 
terizes the inversion asymmetry of the confining potential 
H HI H [13, H, [2g[. This reflects the fact that the 
coefficient a in Eq. (|12p is essentially independent of the 
geometry of the quasi-2D electron systems. In hole sys- 
tems the £ z dependence of spin splitting can be reversed, 
i.e., a large field £ z can give rise to a small spin splitting 
and vice versa [30]. This is due to the fact that the HH- 
LH splittings A hl in Eq. (|15|) depend on the geometry of 
the quasi-2D system that can be tuned by means of £ z . 
Indeed, the implicit dependence of /3{ l on £ z can be such 
that not only, it cancels the explicit £ z dependence in Eq. 
(I14[) . but it can even result in an inverse dependence of 
spin splitting on £ z . 

Experimentally, the spin splitting gives rise to a dif- 
ference Ap = p + — p_ between the densities p+ and 
p_ in the two spin subbands Q. As an example, we 
show in Fig.[2]the measured and calculated spin splitting 
Ap/p for a quasi-2D hole system in a Be-doped (001) 
GaAs-Alo.3Gao.7As single hetero junction at constant to- 
tal density p = p + + p_ = 1.84 x 10 11 cm -2 , where 
the asymmetry was tuned by applying an electric field 
£ z perpendicular to the quasi-2D system using front and 
back gates [30j . It can be seen that the spin splitting 
is reduced when £ z is increased which is opposite to the 
behavior observed in electron systems [24| . 



V. MULTIPOLE EXPANSION OF THE SPIN 
DENSITY MATRIX 

We can gain a more systematic understanding of the 
unusual spin-dependent phenomena in hole systems from 
the spin density matrix which is the fundamental object 
providing a complete description of the system. Neglect- 
ing the orbital degrees of freedom, spin-1/2 electron sys- 
tems are characterized by a 2 x 2 spin density matrix 
p, whereas in spin-3/2 hole systems p becomes a 4 x 4 
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FIG. 4: Spin splitting Ap/p for a quasi-2D hole system in 
a (001) GaAs-Alo.3Gao.7As single hetero junction at constant 
total density p = 1.84 x 10 11 cm -2 , where the asymmetry 
was tuned by applying an electric field £± perpendicular to 
the quasi-2D system [301 ]. The inset shows the calculated HH- 
LH splitting Ai{ that decreases as a function of S z . 

matrix. The dominant character of the occupied eigen- 
states of hole systems depends on the quantization axis 
of the underlying j = 3/2 basis functions. Obviously, 
observable quantities such as the spin polarization may 
not depend on this choice. Therefore, it is necessary to 
formulate the spin density matrix in a way such that ob- 
servable quantities can be calculated independent of the 
particular choice for the basis functions that are used. 
Using the theory of invariants [J, [3l| it is possible to de- 
rive an invariant decomposition of the spin density ma- 
trix of j = 3/2 hole systems that can be interpreted as a 
multipole expansion [321 ] . These multipoles indeed have 
the desired property that they can be evaluated and in- 
terpreted independent of the particular choice of basis 
functions. 

Quite generally, neglecting small terms with cubic sym- 
metry, the spin density matrix p for systems with spin j 
can be decomposed as follows [32| 

p = Y,Pi M i ( 16a ) 

j=0 

= po M + S ■ Mi + Q ■ M 2 + O ■ M 3 + . . . . (16b) 

The 2j + 1 quantities Mj , which we refer to as multipoles, 
are spherical tensors that transform according to the ir- 
reducible representations Vj of the point group SU(2). 
They are analogs of the vector of Pauli spin matrices 
familiar in electron systems. Each Mj is a (2j + 1)- 
dimensional vector with components Mj m , — j < m < j, 
which are (2j + 1) x (2j + 1) matrices. These matrices are 
orthonormal in the sense that tr(Mj m My m r) = 5yy5 mm r. 
They can be calculated using standard angular momen- 
tum theory [33[. M )m for spin- 1/2 and spin-3/2 systems 
have been tabulated in Ref. [12] • The multipole moments 
Pj give the weight of each multipole in the density ma- 
trix. Each p. is a (2j + l)-dimensional vector with scalar 



components p )m . The components pj m of the moments Pj 
can be obtained from 

P* m = tY(M jmP ), (17) 

and the magnitude of each moment is given by 

m 

In Eq. ([T6|l one can thus interpret the vectors Mj as 
multipole operators that "measure" the moments Pj of 
the system. For simplicity we assume that all orbital 
degrees of freedom are integrated over, so that p depends 
only on the spin indices and does not depend on the 
wave vector k. Nevertheless, the coupling of spin and 
k is present in the Hamiltonian through the spin-orbit 
interaction. We note that the dot product for spherical 
tensors appearing in Eq. (|16[) is defined in Ref. [331 ] . 

The moments Pj provide a set of independent param- 
eters characterizing the matrix p (32j. The monopole 
po is identified with the carrier density, while the dipole 
p 1 = S corresponds to the spin polarization or Bloch 
vector at B > 0. These are the only moments present in 
the density matrix of conduction electrons. In j = 3/2 
hole systems, a quadrupole p 2 = Q is also present, which 
reflects the splitting between the HH and LH states dis- 
cussed in the preceding sections. The octupole p 3 = O 
is a unique feature of j = 3/2 systems at B > 0. 

As discussed in Sec. [HI the dynamics of j = 3/2 hole 
systems is characterized by the 4x4 Luttinger Hamilto- 
nian H.L- Obviously, the spherical tensor operators dis- 
cussed here provide the most natural language to formu- 
late the spherical approximation 7i s of Hl, i.e., instead 
of Eq fit w e may write Tt s analogous to Eq. (fT!)]) in the 
form [a. l34| 

3 3 

H s = =^ aj /Cj -Mj (19a) 

j=0 j=0 

= K -M -2V5khbKi -Mi 

+V6^K 2 -M 2 + (K 3 -M 3 + H C . (19b) 
too 

The tensor operators TCj (like the vectors Mj) are tabu- 
lated in Ref. (3^]. The first term in Eq. (pb|) . which we 
interpret as a monopole, is equal to the first term in Eq. 
(J4j) . The second term in Eq. (|19b|) . which corresponds 
to a dipole, is the (spherically symmetric) Zeeman term 
(|8l). The third term in Eq. (j!9b[) . which corresponds to a 
quadrupole, is equal to the second term in This term 
is responsible for the HH-LH splitting. The fourth term 
in Eq. (j!9b[) corresponds to an octupole, which depends 
on k and B. The prefactor £ of this term is typ ically very 
small so that normally it can be neglected [32t[34|. Obvi- 
ously Tii = H.3 = for zero magnetic field. We note that 
in Eq. ^ we considered only the terms with the lowest 
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order of k and B so that we have /Co = k 2 and TCi = B. 
Higher-order terms as well as, e.g., strain-induced terms 
can be classified in the same way. For example, the most 
important strain-induced term is a quadrupole field (in- 
dependent of k [32|). Finally, the cubic term H c reads 

n c = - A- V6( l3 - 72) {[/Ca x M 2 ] 4> _ 4 

+ m [jc 2 x M 2 ] i + [K 2 x M 2 ] 44 }, (20) 

where the product of spherical tensors is defined in 
Ref. [33|. Our prefactor of H c differs slightly from the 
one given m Ref. |[ due to the different normalization of 
the tensor operators 7C 2 and M 2 adopted from Ref. [32| . 

VI. MULTIPOLE MOMENTS INDUCED BY A 
MAGNETIC FIELD 

To gain a qualitative understanding of the relevance of 
the spin multipoles , it is helpful to study the response 
of a 2D hole systems to an external magnetic field B. 
For zero magnetic field, the system is characterized by 
a large quadrupole moment close to its maximum value. 
This result reflects the fact that the HH-LH splitting is 
the most important effect in a 2D HH system and the 
multipole expansion (| 1 6|) is the natural language to de- 
scribe this effect. A perpendicular magnetic field results 
in a large dipole moment (i.e., a large spin polarization), 
as expected from Fig. [1] Of course, for a perpendicu- 
lar magnetic field we also need to take into account the 
formation of Landau levels. A low-density 2D electron 
system can become fully spin-polarized without forming 
Landau levels when an in-plane magnetic field is applied 
[TtI . I35I ] . When the Zeeman energy becomes larger then 
the (B dependent) Fermi energy, the minority spin sub- 
band is completely depopulated and the system becomes 
fully spin-polarized. In 2D HH systems, on the other 
hand, the spin polarization is generally suppressed when 
an in-plane magnetic field is applied [32l |36| . As an ex- 
ample, Fig. [5] shows the calculated response of a 2D HH 
system in a symmetric (100) GaAs-Alo.3Gao.7As quan- 
tum well with hole density p=5x 10 10 cm^ 2 and well 
width w = 150 A. 

We can see in Fig. [5] that the complete depopulation of 
the HH minority spin subband due to B\\ does not imply 
full spin polarization of the system [3^, [Hj]. In Fig. [5J 
the minority spin subband is completely depopulated at 
B D w 23.9 T while S{B D ) « 0.144. Here we use the 
tilde to indicate that we have normalized p\ = S with 
respect to the total 2D density 2po- We see that the spin 
polarization <S(Sm) is always much smaller than 3/2V5 ~ 

0.67, the value of S in a fully spin-polarized HH system 
and it is unrelated with the depopulation of the minority 
spin subband [Fig. Ufa)]. Most surprisingly we even have 
a sign reversal of the spin polarization vector at B» sa 
18.4 T which is a unique feature of 2D HH systems [36] ]. 
The derivative of S(Bu) is discontinuous at B» = Bq- 
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FIG. 5: (a) Spin subband density p- of the HH minority 
spin subband and (b-d) normalized multipole moments pi as 
a function of the in-plane magnetic field B\\ calculated self- 
consistently for a symmetric (100) GaAs-Alo.3Gao.7As quan- 
tum well with hole density p = 5 x 10 10 cm -2 and well width 
w — 150 A. Adapted from Ref. H3). 



At Si I = the quadrupole moment Q is slightly smaller 
than 1/2, the value of Q in a pure HH system. This is 
a consequence of the ky -induced HH-LH mixing which 
was fully taken into account in Fig. O For By > we 

observe only a small decrease of Q. This is due to the 
fact that the HH-LH splitting E[ l - E{ » 6.7 meV [i.e., 
the denominator in Eq. ([9])] is the largest energy scale in 
the system so that the HH states have a frozen angular 
momentum perpendicular to the 2D plane, as illustrated 
in Fig. [TJ For comparison, we note that the Zeeman 
energy splitting of the HH subband at Bp is ~ 0.4 meV. 

It is remarkable that the octupole moment O at B\\ = 
Bo is close to l/\/2 w 0.71, the largest possible value 
of f>3 in a 2D HH system. This value is essentially in- 
dependent of whether we use the Luttinger Hamiltonian 
(|19[) without or with the octupole term proportional to £. 
This result reflects the fact that, unlike the simpler case 
of spin-1/2 electron systems, we cannot establish a sim- 
ple one-to-one correspondence between the multipoles in 
the Hamiltonian (fill)) and the multipoles of the same de- 
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gree in the spin density matrix (JT6J) . On the other hand, 
these findings suggest that an in-plane magnetic field B» 
provides an efficient tool to study 2D HH systems with 
a large octupole moment but with a small dipole mo- 
ment (i.e., with a small spin polarization). Note that we 
can use an in-plane magnetic field to obtain a significant 
out-of-plane spin polarization using the anisotropic Zee- 
man splitting (the component g* z of the g factor) on low- 
symmetry surfaces such as the (113) surface, see Fig. [2] 
(Ref. [H|). 



VII. SPIN PRECESSION OF HOLES 

A. General Analysis 

Before discussing hole spin precession, we will review 
briefly Larmor precession in electron systems [33] • It 
is well known that the dipole moment (spin polariza- 
tion) of a spin-1/2 system performs a precessing motion 
in an external magnetic field B. It can be derived via 
the Heisenberg equation of motion (HEM) for the spin 
operator = cr propagating due to a Zeeman term 



Hi = 29* ^ B <T 



B. 



dcr 

~dt 



ih 



[a, Hi] = g* pb cr x B . 



(21) 



where cr at time t = is the vector of Pauli spin matri- 
ces. If the only spin-dependent term of the full Hamilto- 
nian is Hi, the HEM (|21[) is valid also for particles with 
j > 1/2 (with cr replaced by the appropriate angular 
momentum matrices for j > 1/2). Taking the expecta- 
tion value of Eq. (f2"Tj) with respect to a state \ip) yields 
(&) = g* fiB (<r) x B, which can be interpreted as Ehren- 
fest's theorem [38j applied to a spin-j system. 

In the context of conventional spin precession [37j it 
appears natural that the right-hand side of Eq. (f2"Tj) can 
be expressed as a linear combination of spin operators <r m . 
This implies that the HEM for the components er m of cr 
are closed. Using the language of multipole moments, we 
see here that for spin-1/2 systems the HEM of the dipole 
Mi = cr is decoupled from the HEM of the monopole 
M Q . Obviously, the HEM for M is trivial, M Q = 0, 
which reflects the conservation of the probability density. 
We also note that Eq. (|2"Tj) preserves the length of the 
Bloch vector (cr), i.e., d\ (cr)\/dt = 0. This is equivalent to 
the statement that the magnitude of the dipole moment 
Pi= S does not depend on time, S(t) = constant, which 
reflects the conservation of energy in the system. (Here 
we ignore scattering and spin relaxation [39j.) 

We proceed to study the corresponding equations of 
motion of the multipoles Mj of spin-3/2 systems. While 
the dipole term in the multipole expansion of the Hamil- 
tonian can always be interpreted as a Zeeman-like term 
with an external or effective magnetic field, no such inter- 
pretation is possible for the higher multipoles in Eq. (fTS)) 
including the quadrupole term. Therefore, the simple 
picture of a spin precessing around an effective Zeeman 



field is not applicable to holes. Nevertheless, the spin 
dynamics of hole systems can be viewed as a precession, 
if precession is understood as a nontrivial periodic mo- 
tion in spin space described by an equation of the type 
= j-[H,S] for a suitably generalized spin operator 
S and spin Hamiltonian H. However, the HEM for the 
different Mj cannot be decoupled which has important 
consequences for the spin precession of hole systems [40j . 
Equation (|19a|) suggests that we study first the HEM 



dM, 1 _ „, n 



(22) 



Making use of Eq. (|19ap . this equation can be decom- 
posed into 



(23) 



where the dot product is between Mj' and fCy. The 
commutator is an antisymmetric tensor product that may 
be further decomposed into multipoles using standard 
angular momentum algebra [HI]: 



[Mj 



Here, 



3,m 



m 

j'm' 



l-(-l) 



O+i'+J) 



(24) 



the C- )mym , 
which the phase convention of Ref. [3 



are Clebsch-Gordan coefficients for 
has been adopted. 
Angular momentum conservation constrains the sum to 
|j - j'l < 3 < j + )' and < 971 < 3- The multipole com- 
ponents Mjot, which appear on the RHS of Eq. are 
those satisfying the condition j + + 3 = odd. The invari- 
ant decomposition of the RHS of Eq. (f2"2"|) is summarized 
in Table HI [41j. As discussed above, a dipole Mi evolving 
under the action of Tii describes the well-known Larmor 



precession of holes [37[. Table U shows, as expected, that 
in this case the RHS of Eq. (|22"|) yields only Mi, i.e., 
the HEM for a dipole, which couples to a magnetic field 
through a Zeeman term, is closed. However, the most re- 
markable entry in the table is the one for M2 propagating 
in time due to a quadrupole H.2, the spin-orbit interaction 
that gives rise to the HH-LH coupling. We see here that 
the HEM for the components of M2 are not closed. A 
quadrupole M2 precessing in a quadrupole field "decays" 
into a dipole Mi and an octupole M3. This implies that 
spin precession of an initially unpolarized system can give 
rise to spin polarization, even though B = 0. [As men- 
tioned above, we use the term spin precession for any 
HEM ([221) with i.i' > 0.] 

We want to determine and interpret the explicit time 
evolution of p. This calculation is most easily carried out 
in the Schrodingcr picture, which reflects the equivalence 
of the Heisenberg and Schrodinger pictures for this prob- 
lem. Our final results below are independent of represen- 
tation. In the absence of external fields and disorder, the 
density matrix satisfies the quantum Liouville equation 



dp 
dt 



{ Pl H}. 



(25) 
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TABLE I: Irreducible representations T>j of SU(2) of the (lin- 
ear combinations of) multipole Mj contained in an invariant 
decomposition of (1/ifr) pMj,7tj/]. 





"Ho 


Hi 


H 2 


m 


M 














Mi 





Pi 






M 2 





T) 2 




©2 


M 3 





v :i 


v 2 





Note the sign difference between the Liouville and 
Heisenberg equations. The formal solution is p(t) = 
e -iHt/hp^ e iHt/h^ wnere e iHt/h j g t - me evom tion op- 
erator (which can often be evaluated in closed form). 

B. Precession of a Single Spin 

We consider first an example where we assume the hole 
spin to be oriented initially along the z-direction, rrij = 
+3/2, so that 

P(t = 0) = ^M 0j0 + ^=M lt0 + ^M 2 , + ^=M 3)0 , (26) 

where the tilde indicates that p has been normalized with 
respect to the total density 2p . This equation demon- 
strates that, in general, the density matrix of holes can- 
not be written simply as the sum of a monopole and a 
dipole. The higher multipoles will be present, too [2ll |. 

We want to restrict ourselves to TL = Ho + H2, i.e., 
B = 0. Tablc|T]shows that p evolves as a combination of a 
dipole, a quadrupole, and an octupole. The implications 
of this fact are best seen by following the motion of the 
Bloch vector S(t) = V5S(t) 

S(t) = § {S [cos 2 (ujt)+c 4 sin 2 (ut)] 
+ ksc(l + c 2 ) sin 2 (wi) 
+ (S x k)2sc sin(wi)cos(u^)}, (27) 

where the hat denotes unit vectors; we have So = S(£ = 
0); c — So • k is the cosine of the angle between So and k; 
s is the sine of the same angle; and n = k — (So ■ k)So is 
the vector orthogonal to So in the (So,k) plane. We 
see that the trajectory S(t) in spin space is indepen- 
dent of the Luttinger parameters; only the frequency 
uj = (Eh — Ei)/2h ~ 7?ifc 2 /mo depends on 7 and k. When 
S is parallel to k (i.e., c = 1), we get S(t) = S , which 
is due to the fact that the initial state is an eigenstate of 
the Hamiltonian. In general, neither the magnitude nor 
the orientation of the Bloch vector are conserved. This is 
illustrated in Fig. [6] showing S(t) for an angle of 60° be- 
tween S and k. Hclicity S • k is conserved, a well-known 
fact about this model, which sheds additional light on 
spin precession in hole systems. Since ^(S • k) = and 
the wave vector is not changing, • k = 0. Therefore, 



z 




FIG. 6: Precession of the Bloch vector S (bold arrow) about 
the wave vector k when the angle between the two is initially 
60°. We have assumed S(t — 0) is parallel to the z axis, and 
k lies in the xz plane. Taken from Ref. (4(j. 

whenever the magnitude of the spin changes, the angle 
between spin and wave vector must change in order to 
preserve the projection of S onto k. As a consequence, no 
nontrivial spin precession occurs when S 1 k (or S || k). 
We note that energy is also conserved for spin precession 
in hole systems. Yet for B > 0, when TL = H0 + TL1+H2, 
it can be shown that, in the general case, energy is trans- 
ferred back and forth between Hi and TL2 as time pro- 
gresses. For an LH spin (rrij = +1/2) we obtain similar 
to Eq. d2Tj) 

S(t) =S [| + § s 2 (2 - 3s 2 ) sin 2 (urf)] 

+ § ksc(l - 3c 2 ) sin 2 (ujt) 

— 3 (So x k) sc sm(wi) cos(uit). (28) 

The analysis above highlights a major advantage of our 
method. Although most of the results presented here can 
be derived also using wave functions, the decomposition 
of p and TL into multipoles makes the symmetry of the 
problem transparent [8j. In particular, the interdepen- 
dence of the multipoles would not be evident if one used 
wave functions. Moreover, the multipoles do not rely on 
a particular choice of basis functions [13] — in that sense 
the multipoles are gauge-invariant. 

C. Holes Turning a Corner at B = 0: Alternating 
Spin Polarization 

We would like to relate our work to recent experiments 
by Grayson et al. [13, HH which demonstrated that a 
high-quality bent heterostructure can be grown on top 
of a pre-cleaved corner substrate that allows one to drive 
the charge carriers around an atomically sharp 90° cor- 
ner. Grayson's experiments were performed on a 2D elec- 
tron system in a GaAs/AlGaAs quantum well. We show 
here that a similar system containing holes gives rise to 
fascinating new physics [40(. The setup is sketched in 
Fig. [7Ja) . We assume B — so that the Hamiltonian is 
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FIG. 7: A bent structure allows holes to be driven around a 
corner [U HH). (a) An idealized representation of the corner 
device, (b) A more realistic model of the corner, as discussed 
in the text. In regions C and L2 the spin polarization oscillates 
as a function of time. Taken from Ref. 401. 



H = Ho + Ti-2 ■ An unpolarizcd HH wave packet travels 
in the 2D channel Li in the +x direction. For simplicity 
we assume that all orbital degrees of freedom are inte- 
grated over, so that p depends only on the spin indices 
and does not depend on the wave vector k. Nevertheless, 
the coupling of spin and k is present in the Hamiltonian 
through the spin-orbit interaction. The magnitudes of 
the normalized moments at t < (i.e., before reaching 
the corner) are S = O = and Q = 1/2. The spin 
quantization axis of the HH states in Li is parallel to the 
z-direction. After the wave packet has passed the corner, 
the HH states are not eigenstates of H2 , their spin quan- 
tization axis being perpendicular to the spin quantization 
axis supported by the quadrupole field in L2. Therefore, 
the quadrupole and octupole moments in L2 oscillate in 
time 



\Q(t)\ 2 = (1/16) + (3/16) cos 2 (w 2 t) 
\0(t)\ 2 = (3/16) sin 2 (Lu z t), 



(29a) 
(29b) 



with precession frequency lu z = (Eh — Ei)/2H ~ 
2jhn 2 /mQW 2 . This frequency can be tailored by vary- 
ing the width w of the 2D channel. For GaAs in the 
spherical approximation we have 7 = 2.58 and we take 
w = 10 nm, yielding a precession period 2tt /uj z 11 ps. 

For the simplified geometry of Fig. UJa) the spin po- 
larization S in L2 remains zero, as required by the con- 
servation of helicity discussed above. However, assuming 
a sharp 90° corner for an ideal 2D system is certainly an 
oversimplification, even when the corresponding quasi 2D 
system has atomically sharp interfaces [42[ . A more real- 
istic treatment can be obtained by modeling the transi- 
tion region between the channels Li and L2 as a sequence 
of two 45° corners as sketched in Fig. UJb). Once again, 
an unpolarized HH wave packet travels in channel Li in 
the +x direction, the initial conditions being the same as 
in the previous example. If the HH wave packet enters 
the central region C at t — 0, we obtain for the squared 
normalized moments in this region 

\S{t)\ 2 = (9/80) sin 2 {uj z t) (30a) 
\Q{t)\ 2 = 1/64+ (15/64) cos 2 (w z i) (30b) 
|e>(^)| 2 = (39/320) sin 2 (w z t). (30c) 



Equation ()30a|l shows that the initially unpolarized hole 
current acquires an alternating spin polarization S(t) 
due to spin precession at B = 0. In Cartesian co- 
ordinates, the Bloch vector in region C reads S(t) = 
[0, — -j^g sin(ui z t), 0]. When the HH wave packet enters 
the channel L2 it continues to precess. We get for the 
spin polarization 

~ —3 2 ( LJ Z T\ , 2 / Wji 

S y = ^ ,_ cos I - j sm(uj z t) + 2 sin(Lu z T) sin ^— — 

Here T is the time required to traverse C which depends 
on the length of C and the magnitude of the in-plane wave 
vector. It determines the fraction of the spin polarization 
in L2 that will be oscillating. If we take the length of C 
to be of the order of the channel width, namely w = 
10 nm, and the initial wave vector kp = 0.1 nm" 1 , we 
get T ~ 0.5 ps. The amplitude of S y in this case is 
approximately 0.1. We omit here the qualitatively similar 
but more complicated expressions for Q(t) and 0(t). We 
note also that the approach in Fig. [7(b) can be further 
extended in a transfer-matrix-like approach in order to 
describe more complicated geometries. 



VIII. SPIN RELAXATION IN HOLE SYSTEMS 



Spin relaxation in spin- 1/2 electron systems has re- 





ceived considerable attention 

l66l |67| . For electrons the spin-orbit interaction 



— 
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can always be represented by a Zeeman-like term f| L3[) 
with an effective wave vector-dependent magnetic field 
B(k) = ftf2(k). The electron spin precesses about this 
field with frequency = |fi(k)|. An important parame- 
ter is the product of the frequency f2 times the momen- 
tum relaxation time r v . In the ballistic (clean) regime no 
scattering occurs, so that Qt p — > 00. The weak scattering 
regime is characterized by fast spin precession and little 
momentum scattering, yielding flr p 3> 1. In the strong 
momentum scattering regime VLt p -C 1. Electron sys- 
tems are often in the strong scattering regime and most 
past work has concentrated on this case. 

As discussed in Sec. [V] for spin-3/2 holes the spin- 
orbit interaction cannot be written as an effective field, 
and spin precession is qualitatively different (Sec. IVII[) . 
Since spin-orbit coupling is more important in the va- 
lence band, hole spin information is lost faster, and the 
relative strengths of spin-orbit coupling and momentum 
scattering can vary. Yet spin relaxation of spin-3/2 holes 
has also been studied to a lesser extent, both experimen- 
tally [H and theoretically (HE M, M, [ED, M, lz| ■ 

Recently we were able to derive a general unifying 
quantitative theory for the return to equilibrium of excess 
spin polarizations in the conduction and valence bands 
of semiconductors brought about by the interplay of spin 
precession and momentum scattering [39| . Spin polariza- 
tion decay in different regimes of momentum scattering 
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in spin-1/2 electron and spin-3/2 hole systems contains 
considerable rich and novel physics. For example, spin 
polarization decay has often been assumed to be propor- 
tional to e~ f / Ts , where t s is referred to as the spin relax- 
ation time. However, if the magnitude of the spin-orbit 
interaction is anisotropic (as is usually the case in systems 
studied experimentally), spin-polarization decay can oc- 
cur even in the absence of momentum scattering. This 
process is characterized by a non-exponential decay and 
is sensitive to the initial conditions, and cannot therefore 
be described by a spin relaxation time. Weak momentum 
scattering introduces a spin relaxation time t s cx t p (un- 
like str ong momentum scattering, which gives the well- 
known [USB] trend t s cx t^ 1 ), yet even in the presence 
of weak momentum scattering a fraction of the polariza- 
tion may survive at long times. Indeed, in the ballistic 
and weak momentum scattering regimes, the concept of 
a spin relaxation time is of very limited applicability and 
in general does not provide an accurate description of the 
physics of spin polarization decay. 

A. Time Evolution of the Density Matrix 

We assume a nonequilibrium spin polarization has 
been generated in a homogeneous, unstructured system 
and study its time evolution in the absence of exter- 
nal fields. The system is described by a density matrix, 
which in principle has matrix elements diagonal and off- 
diagonal in momentum space. Since the spin operator 
is diagonal in the wave vector k, we will only be con- 
cerned with the part of the density matrix diagonal in 
momentum space, which is denoted by pk- 

The spin density is given by (S) = tr(Spk) = tr(Spfc), 
where S is the spin operator, the trace includes k and 
spin, and the overline represents averaging over directions 
in momentum space. Only the isotropic part pk of the 
density matrix is responsible for spin population decay 
(46| . It is therefore convenient to divide pk into pk = 
Pk + .9k, where <?k is the anisotropic part of pk. Based on 
the quantum Liouville equation, we obtain an equation 
describing the time evolution of pk (Ref. 56}), which in 
turn is split into a set of equations for pk and gk similar 
to those found by Pikus and Titkov [4a ]: 

^ + ^I^kT-0, (32a) 
dgk i r _ , ffk dp k i _ 

^T + H [n ' 9k] + V p = -dt-h [n > Pk] - (32b) 

Here, H is, in general, the full Hamiltonian; yet the 
spin-diagonal part of 7i commutes with pk so that it 
is not relevant for spin relaxation. We assume elastic, 
spin-independent scattering, implying that the collision 
term involving pk vanishes [46]. For isotropic (s wave) 
scattering due to, e.g., screened impurities the remain- 
der is proportional to the inverse of the scalar momen- 
tum relaxation time l/r p (Ref. [zlj]). A solution to Eq. 



(|32bp can be obtained by making the transformation 
9k = e~ mt l n g^ e l ' Ht / h , which is analogous to the cus- 
tomary switch to the Heisenberg picture. Substituting 
this solution into Eq. (|52"a|) yields [39| 

^Pt + J_ f e -(t-*Vr, e -iH{t-t')/H [ Wj p k i t ')] e tH(t-t')/h 
at hT p J 

= - ~ e- t/T ? e - mt / h [W,pk ] e mt / h , 

where pko is the initial value pk(£ = 0). This equation de- 
scribes the precession-induced decay of spin polarization 
in all regimes of momentum scattering for both spin-1/2 
electrons and spin-3/2 holes in semiconductors. It does 
not anticipate any particular form of spin polarization 
decay, such as exponential decay. 

The form of the initial density matrix pko is impor- 
tant and lies at the root of the physics discussed in the 
remainder of this section. In general pko has two contri- 
butions, pko = Pko + Pko- The component p k0 commutes 
with H. and pjjr is simply the remainder. pJL is a matrix 
that is parallel to the Hamiltonian, and represents the 
fraction of the initial spin polarization that does not pro- 
cess, or alternatively the fraction of the initial spins that 
are in eigenstates of the Hamiltonian. p^ is orthogonal 
to the Hamiltonian, and represents the fraction of the 
initial spin polarization that does precess. 

B. Regimes of Spin Polarization Decay in Spin-3/2 
Hole Systems 

In general, Eq. (|3"3"]) allows one to distinguish and de- 
scribe several different regimes of spin polarization decay 
by comparing the precession frequency Q (k) with the mo- 
mentum relaxation time r p [3^ ] . In the following we con- 
sider spin-3 /2 holes described by the Luttinger Hamilto- 
nian (TTJ. We work first in the spherical approximation, 
Eq. (HU). 

1. Exponential Decay in the Strong Momentum Scattering 
Regime 

A solution to Eq. (|33[) characterizing relaxation is un- 
derstood as exponential decay of the form 

Pk(t)^e~ F ^p k0 , (34) 

where T s is generally a second-rank tensor that repre- 
sents the inverse of the spin relaxation time r s . Such a 
simple solution of Eq. (|33[) does not exist in general, but 
for strong momentum scattering (£lr p <C 1) the RHS of 
Eq. (j3U|) can be neglected. Then substituting for pk and 
TL in Eq. (|33[) yields an exponential decay of the spin 
population with a relaxation time t s cx t^ 1 . This trend 
is well-known for Dyakonov-Perel spin relaxation in elec- 
tron systems [3 . 
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For spin-3/2 holes we get T s = r^ 1 1, showing that 
(for a given wave vector) the relaxation times for all spin 
components are equal, 

± = itfT p = U^) 2 r p , (35) 

where O(k) = (Elh — £hh)/ ft = 2h^k 2 /m,Q. The spin re- 
laxation times r s for HH states become effectively smaller 
than t s for LH states if we take into account that HH 
states are normally characterized by a larger Fermi wave 
vector k F . Despite the qualitatively different spin pre- 
cession, the situation is overall rather similar to electron 
spin relaxation and can be explained in terms of the same 
random walk picture familiar from the study of electron 
spin relaxation [3, EB| • 
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FIG. 8: Incomplete spin dephasing of (a) electron spins in 
the fc 3 -Dresselhaus model and (b) of heavy-hole spins (solid 
line) and light-hole spins (dashed line) in bulk GaAs in the 
ballistic limit. The vertical axis shows the normalized spin 
polarization S(t) ■ So/|So|. The initial spin polarization So 
is assumed to point along [001]. In (a) the Fermi energy is 
Ef = 5.7 meV; in (b) we used Ef = —1.0 meV. Note the 
different time scales in (a) and (b). 



2. Ballistic Regime in the Spherical Approximation 

In the opposite limit t p — > oo the second term on the 
LHS of Eq. d33j can be neglected. Then Eq. (J33J) is solved 
by 



Pk(t)=plo + ^ mt/h PL^ Ht/h (36) 

which describes a spin precession of the initial spin polar- 
ization pko- In the spherical approximation (fT9f for the 
Luttingcr Hamiltonian, an initial spin polarization will 
oscillate indefinitely since f2 is the same for all holes on 
the Fermi surface. 



3. Weak Momentum Scattering Regime in the Spherical 
Approximation 

In the regime of weak momentum scattering the solu- 
tion to Eq. (|33|) may be written approximately as 

Pk(t) = p1 + e- tlr * e-w/Kp^eWI*, (37) 

Since the momentum scattering rate 1/t p is small, the 
term under the over line is taken to lowest order in l/r p . 
The second term on the RHS of Eq. ([BT]) describes 
damped oscillations with amplitude decaying exponen- 
tially on a scale cx t p . The fraction S k0 of the spin polar- 
iazarion corresponding to p k0 survives at long times. We 

can determine S k0 by averaging Eqs. (|27|) and (|28f over 

I 6 

time t and directions of k, showing that S k0 = ^jSko 
for both HH and LH states. This result does not depend 
on the Luttinger parameters or the Fermi wave vector 
and will therefore be the same in any system described 
by the Luttinger Hamiltonian ri s . The remaining polar- 
ization S yn d ecays via spin-flip scattering as discussed in 

Refs. (zalzl. 



4- Cubic- Symmetry Terms and Dephasing 

Dephasing is introduced if the term ri c with cubic sym- 
metry is included in the Luttinger Hamiltonian. The 
cubic-symmetry terms contained in Eq. (|20|) are usually 
neglected in charge and spin transport without a signif- 
icant loss of accuracy. Due to the presence of H. c , the 
energy dispersion relations and therefore |fi(k)| depend 
on the direction of k. Spins on the Fermi surface thus 
precess with incommensurjable frequencies and once they 
are out of phase they never all get in phase again. Even in 
the ballistic limit this process results in a non-exponential 
spin decay [H, [75| with a characteristic time T<j cx 
referred to as the dephasing time Td- In general (in par- 
ticular for 3D systems), an intermediate situation is re- 
alized where the spin polarization is reduced because of 
dephasing, but it remains finite. The surviving part is 
identified with /? k0 in the initial density matrix. This 
process is referred to as incomplete spin dephasing. 

Our numerical calculations exemplified in Fig. [5] show 
the incomplete dephasing of electrons and holes in bulk 
GaAs. For electrons, dephasing is caused by the k 3 - 
Dresselhaus model [76|. At long times the initial spin 
polarization settles to a value w 0.33, which is indepen- 
dent of any system parameters, including the spin-orbit 
constant. For holes, the initial spin polarization falls to a 
fraction much higher than in the electron case. It decays 
more slowly for the LHs, for which the Fermi surface is 
nearly spherical, than for the HHs, for which the Fermi 
surface deviates significantly from a sphere. Note also 
that for a given Fermi energy the HH Fermi wave vector 
k™ is much larger than k F K so that the HH states pre- 
cess faster than the LH states. At long times the spin 
polarization settles to a value w 0.59 for the HH states 
and ~ 0.70 for the LH states. These nonuniversal val- 
ues differ slightly from the universal value 6/10 of the 
surviving spin polarization obtained in the spherical ap- 
proximation. 

In Fig. [5] we assumed that the initial spin polarization 
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is isotropic in k space. It is known [73] that optically 
excited spin distributions are highly anisotropic. Indeed, 
it turns out that optically oriented heavy or light holes 
in 3D do not precess at all (i.e., p^ = 0) which is similar 
to 2D electrons in a symmetric QW on a [110] surface 

mm. 

IX. CONCLUSIONS 

We have reviewed spin orientation in semiconductor 
hole systems that are characterized by an effective spin 
j = 3/2. We showed that the Zeeman splitting and 
Rashba spin splitting in hole systems are qualitatively 
different from their counterparts in electron systems. A 
systematic understanding of the unusual spin-dependent 
phenomena in hole systems can be gained using a multi- 
pole expansion of the spin density matrix. As an example 



we discussed the spin precession in hole systems that can 
give rise to an alternating spin polarization. Finally we 
discussed the qualitatively different regimes of hole spin 
polarization decay in clean and dirty samples. 
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